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Normal modes analysis of the microscopic dynamics in hard discs 
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We estimate numerically the normal modes of the free energy in a glass of hard discs. We 
observe that, near the glass transition or after a rapid quench deep in the glass phase, the density 
of states (i) is characteristic of a marginally stable structure, in particular it displays a frequency 
scale io* ~ p^ , where p is the pressure and (ii) gives a faithful representation of the short-time 
dynamics. This brings further evidences that the boson peak near the glass transition corresponds 
to the relaxation of marginal modes of a weakly-coordinated structure, and implies that the mean 



square displacement in the glass phase is anomalously large and goes as (5R^) 
that we check numerically. 

PACS numbers: 61.43.-j, 64.70.Pf, 67.55.Jd 
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There is still no accepted description of why the dy- 
namics of liquids rapidly slows down as the temperature 
is lowered [l] . Goldstein [^ proposed that the glass tran- 
sition is related to a change in the topology of the free 
energy landscape. In this view, below some temperature 
meta-stable states appear, the dynamics becomes acti- 
vated and rapidly slows down. A similar scenario occurs 
in mean field glass models [3, |j, [^] , which share formal 
similarity with the Mode Coupling Theory of liquids l(| . 
The presence of an elastic instability, corresponding to 
the onset of activation, is also predicted in systems were 
particles are deposited randomly (corresponding to an 
infinite temperature) , as the density is changed [7| . Nev- 
ertheless, the real-space interpretation of those views is 
not clear. Empirically, the presence of an elastic instabil- 
ity near the glass transition is suggested by Raman and 
neutron spectra, which display a peak, the "Boson Peak" , 
typically shifting toward zero-frequency as the glass is 
heated toward the glass transition^. This observation 
is also supported by numerical simulations '9J. The na- 
ture of the modes forming the peak and how they depend 
on the structure of the glass are long-standing questions 
[10|, on which recent progress was made in a-thermal as- 
semblies of particles [ll|, [l2| . Elucidating these questions 
at finite temperature, in particular in the vicinity of the 
glass transition, is necessary to make progress in our mi- 
croscopic understanding of the glass transition. Indeed 
these modes, which characterize the short-term dynam- 
ics, presumably affect how glasses flow on much longer 
time-scales, as supported by the observations that the 
short term dynamics (i) correlate with the fragility of 
the glass [13] and (ii) is a good local predictor of the 
propensity (the tendency of a region to flow) [1J|. To 
study these modes, various numerical technics can be 
used, such as the instantaneous normal modes analysis 
[15| or the analysis of the vibrational spectrum of inher- 
ent structures [ly|. Nevertheless, these approaches con- 
sider the energy rather than the free energy landscape. 
Such analysis are not adapted for colloidal systems or 
shaken grains [ITI . Il8j or when non-linearities of the in- 



teracting potential are important, which is presumably 
often the case near the glass transition p^ . 

In this Letter, we study the microscopic dynamics of 
hard discs in the vicinity of the glass transition and deep 
in glass phase, via an analysis of the free energy land- 
scape. To achieve it, we use a recent technic developed 
to infer the normal modes of the free energy [20|. We 
show that, (i) after a rapid quench in the glass phase or 
near the glass transition, the density of normal modes 
is characteristic of a marginally rigid solid, in agreement 
with a former analysis of the microscopic structure of 
the glass [20|, (ii) throughout the spectrum, the com- 
puted frequency of a mode gives the correct estimate for 
the characteristic time and amplitude of its fluctuations. 
This demonstrates that our mode analysis performs well 
in the vicinity of the glass transition, (iii) These results 
justify the presence of a boson peak and an anomalously 
large and slow microscopic relaxation. In particular, the 
mean square displacement must vary as {SR'^) ~ p~^'^, 
in good agreement with our simulations. This signifies 
that, deep in the glass phase, the amplitude of the parti- 
cles trajectories become infinitely larger than the instan- 
taneous size of the cages. 

As the viscosity of a super-cooled liquid increases, the 
dynamics become intermittent [2ll . [23 . [23| : there are long 
quiet periods where particles are rattling around their 
average positions, interrupted by rapid rearrangements 
toward new configurations. In a simulation, this inter- 
mittent dynamics can be followed via the self density 
correlation function C{q,t,t^) = (e»«-(«j(*+*™)-^j(*™)))^-, 
where the average is made on all particles, but not on 
time, and where Rj{t) indicates the position of particle j 
at time t. This observable displays long plateaus, inter- 
rupted by sudden jumps leading to new plateaus [20ll2]|. 
We shall refer to these plateaus as meta-stable states. In 
the reference [21| , the transition between states was stud- 
ied. Here we focus on their inner dynamics. 

As shown in 20], in meta-stable states the Gibbs free 



energy of a hard sphere system can be approximated as: 



g = -hTj2Hh 
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where the sum is made on all pairs of particles ij in con- 
tact, i.e colliding with each other in a specific meta-stable 



state. T is the temperature, {hi 



\Ri — Ri 



{di 



crj)/2 is the average gap between particles measured in a 
given meta-stable state, and Oi is the diameter of parti- 
cle i. Eq.([T|) draws an analogy between the free energy of 
hard spheres and the energy of a system of logarithmic 
springs. Eq.lJT]) can be expanded for small displacements 
around the average particles position bRi = Ri — (Ri). 
This enables to define the dynamical matrix A4 |24| . 
which describes how the average particles positions re- 
spond to a small external force. The eigenvectors \5R'^) 
of A4 are the vibrational modes [25| . For particles of unit 
mass, we define as usual their frequency u! as the square 
root of the corresponding eigenvalue of A4 . In what fol- 
lows u! is used to label the normal modes. We shall see 
below that this label corresponds, within a good approxi- 
mation, to the physical oscillation frequency of the mode. 

In practice, we use an event-driven code to simulate 
a two-dimensional bidisperse hard discs system [201 . |2]| 
and compute Ai . Half of the particles have a unit diam- 
eter, while the others a diameter of 1.4 and their mass 
is m = 1. We choose k^T as our unit of energy. We can 
equilibrate this system up to < (/)o ~ 0.79. To generate 
larger packing fractions and enter in the aging regime, 
we proceed as follows: we decompress jammed configu- 
rations of [2(| of packing fraction (j)c ~ 0.83 by decreasing 
the particle diameters by some fraction, and assign a ran- 
dom velocity to every particle. As the dynamics run, we 
compute the self-density correlation function and deter- 
mine its plateaus or meta-stable states. For each of them, 
we consider a time interval [t,t + ii] corresponding to a 
single state, where ti is chosen to be much larger than 
the average collisional time of a contact Tc. In what fol- 
lows ii corresponds typically to 100 Tc but our results are 
robust to other choices as long as ii 3> Tc. The network 
of contacts and the average particle positions (Ri) are 
computed during this interval. This enables to compute 
g and A4 numerically. 

To compute the normal modes and the density of 
states, we diagonalize M. Results are shown in Fig.([T]). 
Curves are labeled by the average contact force (/) = 
{i/{hij))ij, where the average is made on all contacts. 
(/) diverges near 0c and scales as the pressure. The dens- 
est system we can equilibrate corresponds to (/) ~ 18. 
For larger (/) the system is a glass. To compare vari- 
ous packing fractions, the frequency axis is rescaled by 
the square root of the typical interaction stiffness fc, of 
order fc^/^ ~ (/) in our units. We observe that: (i) for 
all (/), D{io) increases rapidly from zero-frequency to 
reach a maximum at some frequency to* , before decaying 
again. In our finite size system the modes constituting 
the low-frequency part of the spectrum are not plane- 
wave like 28i. This can be inferred from the inset of 
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FIG. 1: Densities of states D'(lo') = {f)D(u) vs. rescaled 
frequency ui' = u!/{f} for different values of (/) in a system 
oi N = 256 particles. We used ti = 5 x 10*. Inset: D{uj')/uj' 
vs. Lo' for {/) = 18. 



Fig.©, where D{ijj) is normalized by its Debye value 
(^ Lj in two dimensions): no plateau can be detected at 
low frequency, where a peak appears at some frequency 
Wbp significantly smaller than lo* . (ii) lo* decreases to- 
ward zero in the rescaled units as (/) diverges. Fig. ([2) 
shows the dependence of lo* with (/). We observe that 
the scaling 
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holds well from the glass transition toward our densest 
packing. Scaling behaviors in the vibrational spectrum, 
together with the fact that the lowest-frequency vibra- 
tions observed are not plane wave like, were first observed 
in elastic non-thermal particles [2y|. This was later ex- 
plained as follows: in a weakly-coordinated network of 
springs at rest, i.e. with a small coordination number z 
close to a critical value Zc — 2d (d is the spatial dimen- 
sion), extended modes, quite different from plane waves, 
appear already in the low-frequency part of the spec- 
trum. The frequency scale of these anomalous modes is 
LO* ~ k^/^{z — Zc) [ll|. In a system of repulsive par- 
ticles with non-zero forces, these modes are shifted to- 
ward a lower frequency [23. To avoid an elastic insta- 
bility, the system must display enough contacts. In the 
marginally stable case where there are just enough con- 
tacts to be rigid, the spectrum displays anomalous modes 
up to zero frequency, but present a characteristic fre- 
quency at LO* ^ if)^^^, such that D{lo*) is of order of its 
typical value A:~^/^. This estimate of D{lo*) agrees with 
Fig-©! where D{lo*) appears to remains of the same or- 
der in rescaled frequency. Thus both the scaling of oo* 
and the observation that anomalous modes are present 
at frequency lo-b-p « lo* support that the glass phase 
lies close to marginal stability, as was already suggested 
by a direct measure of the coordination in the glass phase 



We check that D{lo) reflects well the dynamics. We 
project the dynamics on the normal modes and define 
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FIG. 2: Characteristic frequency lo* as defined in the text vs 

average force (/). Slashed line corresponds to the fit w* ~ 



for each frequency w: 



(3) 



where |(5_R(t+iu,)) is the displacement field from the aver- 
age configuration, (6R°-\6T&) = ^^ &R1 ■ SR^, and where 
the average is made on all time segments [t^ , t^, + i] in- 
cluded in a meta-stable state. If the modes were longitu- 
dinal plane-waves, Ca,(t) would describe the relaxation of 
density fluctuations. Here this is not at all the case. We 
compute the amplitude A(lu) of a mode fluctuation de- 
fined as {A'^{uj)) = Cui{0), and a characteristic time t{uj) 
of the mode relaxation, defined as the earliest time at 
which C^{t) has decayed by 10%: C^{t{uj)) = 0.9C„(0). 
Taking other thresholds (30% or 50%) leads to the same 
scaling result but reduce our statistics. The dependence 
of these quantities with frequency are respectively shown 
in Fig.(l3]-a) and Fig.(l3]-b) for three meta-stable states at 
different pressure. Two configurations are in the glass 
phase and one in the super-cooled liquid phase. In all 
cases, these quantities were computed for each mode of 
the spectrum. We observe that the modes display weakly 
damped oscillations, whose amplitude and period follow: 
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These results hold true even for the low-frequency part 
of the spectrum, although more scattering is found there 
[2 9] . Thus, even for the lowest-frequency part of the spec- 
trum, the non-linearities do not appear to dominate the 
mode dynamics in a dense disordered assembly of col- 
liding hard spheres, although the interacting potential is 
discontinuous. 

This observation implies that -D(w) gives a rather faith- 
ful distribution of relaxation time scales of the micro- 
scopic dynamics. This allows us to identify the peak 
apparent in the inset of Fig.([T]) as the Boson Peak, which 
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FIG. 3: (a) Average squared amplitude of the modes {A (lo)) 
vs LO at various packing fractions in a system of 256 particles, 
both in the glass phase ((/) = 260 and {/) = 50) or in the 
super-cooled liquid ((/) = 18). Each point corresponds to one 
mode. The dashed line corresponds to the fit {A^({jS)) ~ XjuP' . 
(b) Relaxation time t(w) of each mode vs uj for the same 
packing fraction. The slashed line corresponds to the relation 
t{lu) ~ 1/uj. 



appears as a similar h ump in Raman or neutron spectra 
in molecular liquids @,[iO|- Near the glass transition, this 
peak appears at a frequency significantly smaller than 
(jj* . The scaling behavior of D{ijj) supports that these 
modes correspond the anomalous modes characterizing 
weakly-coordinated structures, whose geometry is stud- 
ied in [231 ■ It would be interesting to go further and test if 
the early-/? relaxation [6[, which characterizes the relax- 
ation on time scales longer than the boson peak, is due to 
the weakly-damped relaxation of anomalous modes even 
softer than the peak, or to other non-linear processes not 
captured by the present analysis. Unfortunately this re- 
laxation is difficult to study numerically using Newtonian 
dynamics, and our present analysis should be performed 
in three-dimensions, where numerical and experimental 
data exist, to test this possibility. 

In what follows we rather focus on the amplitude of 
the short-time relaxation, in particular the particle mean- 
square displacements (SRf). Assuming that the dynam- 
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FIG. 4; Mean square displacement (SR^) versus average con- 
tact force (/) for A'' = 1024 (circles) and A'' = 256 (squares) 
particles. Diamonds (red points) correspond to equilibrated 
systems of A'' = 256 particle. Slashed line corresponds to the 
best fit of the form {SR^) ~ (,f)"^/^ 



ics of different modes is independent, one gets: 



(<5i?,')=5](A2(.;))Ji?^ 



(6) 



We then average on all particles and define {SR ) — 

1/A^ ^i{SRi ) where N is the system size. Using the 

modes normalization {5Ri{uj) )i — \/N and applying 
Eqs.® lead to: 
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Since D{uj) reaches the order of its typical value 1/vfc ~ 
l/(/) for cu > uj*, the last integral is dominated by the 
lowest bound and one gets: 
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near maximum packing, the amplitude of particles mo- 
tions becomes infinitely smaller in the crystal, where it 
is simply of the order of the inter-particle gap h. In the 
hard sphere glass, particles diffuse much more than the 
size of their instantaneous cage h, and the size of the re- 
gion they visit in a meta-stable state is at least of order 

To check this prediction, we consider various meta- 
stables states, and measure the mean square displace- 
ment {SR'^) by averaging on some long time ii, see 
Fig. (HI . Our numerical result agrees well with ^. As 
for Eq.([2]), the scaling prediction captures the behavior 
also near the glass transition, where the system can be 
equilibrated. 

We have used and tested a free energy-based approach 
that allows us to perform a mode analysis at non zero- 
temperature, in particular in the vicinity of the glass 
transition. We have observed that, near the glass tran- 
sition or after a rapid quench, the spectrum of the free 
energy is characteristic of a marginally rigid solid, as con- 
firmed further by the scaling of the amplitude of the mean 
square displacement in the glass phase. The observation 
that on the time scales we can probe the system remains 
near the onset of rigidity supports the proposition, ex- 
pressed originally by Goldstein in terms of free energy 
landscape, that the dynamics dramatically slows down 
once meta-stable states appear. Our analysis provides 
a spatial interpretation of this landscape picture, where 
the elastic instability is governed at a microscopic level 
by coordination and contact forces, and where the soft 
modes at play are collective rearrangements, whose char- 
acteristic length scale is limited near the glass transition 
and diverges near the maximum packing of hard spheres 
(ll| . Finally, the observation that the microscopic dy- 
namics, and activation events [2l|, are dominated by soft 
modes is consistent with the observation that the short 
term dynamics is a good indicator of propensity |14| | . It 
suggests that correlations between coordination, contact 
force, and propensity may exist; and that soft modes may 
be the elementary objects to consider to describe micro- 
scopically activation. 



It is easy to show from Eq. ([7|) that the scaling ^ be- 
comes an equality if D(uj) does not increase more rapidly 
than linearly between and lu* [30|- The bound ([8|) 
holds in any dimension d > 2, as long as the marginal 
stability of the structure and the harmonic behavior of 
the modes is assumed. Eq.([5]) must be compared with 
the fluctuations of a crystal: (SR^) ^ h'^ ^ (/)~^ (with 
logA^ corrections in two dimensions). In other words. 
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